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We review some experimental and theoretical results on the metal-to-insulator transition (MIT) 
observed at zero magnetic field (B — 0) in several two-dimensional electron systems (2DES). Scaling 
of the conductance and magnetic field dependence of the conductance provide convincing evidence 
that the MIT is driven by Coulomb interactions among the carriers and is dramatically sensitive to 
spin polarization of the carriers. 



I. INTRODUCTION 

Landauer's early work on scattering of particles and 
waves in random media is still important today, be- 
cause after decades of work, the physics of transport in 
disordered systems is still keeping us entertained with un- 
expected results. The generic "conductivity" that arises 
from ham-handedly averaging microscopic details and 
sweeping Coulomb interactions under the rug of effec- 
tive mass is a poor description of transport in many cir- 
cumstances. In the past two decades there have been 
several discoveries in experiments (much of it spurred or 
explained theoretically by Landauer and his co-workers) 
that have flown dramatically in the face of the classi- 
cal conductivity. Response from quantum mechanically 
phase-coherent carrier excitations have appeared in stud- 
ies of metals and semiconductors in both the relatively 
clean and very dirty limits of impurity scattering. In high 
magnetic fields the re-writing of Ohm's law has been es- 
pecially dramatic. Noise from fluctuations in the trans- 
port parameters have behaved in ways that contradict 
the simple models that lead to the generic resistivity. 
Q Finally, Coulomb interactions are appearing as large 
perturbations of the non-interacting carrier models that 
have dominated thinking for two decades. 

Since the appearance of the scaling theory for non- 
interacting (Fermi-liquid) particles most have ac- 
cepted the result that there is no minimum metallic con- 
ductivity in three dimensions (3D) and no metallic be- 
havior at all in two dimensions (2D). The non- interacting 
models were developed in terms of a scaling function 



that related dimensionless conductance g = G/(e 2 /h) (G 
is the conductance of the sample in l/il) to the size L 
of the (hyper-cube) sample. For very low conductance 
(g <C 1), the transport is exponentially reduced as L in- 
creases because the carriers are localized and must hop or 
tunnel from site to site. For g ^> 1, Ohm's law is assumed 
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FIG. 1. Schematic theoretical plot of f3 with calculated 
results (solid lines) for noninteracting particles in the limits 
/3 3> 1 or j3 <C 1 and plausibility arguments (dotted lines) con- 
necting them. Alternative versions (dot-dash line) for d = 2 
have been proposed for interacting particles. 

to hold: g cx L d ~ 2 for cubes of dimension d. In weak 
disorder, perturbative calculations from non-interacting 
models of carrier transport lead to reduction of (3 indi- 
cating a (perfectly plausible) trend toward lower conduc- 
tance and eventually to an insulator as the amount of 
disorder increased (illustrated in figure Q). According to 
this picture, a 3D system is a conductor if the Fermi en- 
ergy Ep is larger than the characteristic amplitude of the 
disorder potential W. This criterion is equivalent to the 
Ioffe-Regel criterion that Uf I > 1, where Uf is the Fermi 
wave-number and I is the mean free path. In 2D, there 
is no metallic behavior for any value of W ^ 0. That is, 
the proposed 2D scaling function (3 < 0. Both procla- 
mations are reasonable at first glance, but in both cases 
the Coulomb interactions among the carriers have largely 
been swept beneath the rug of effective mass and (of- 
ten ignored) Fermi liquid parameters. Subsequent work, 
however, has brought to question the validity of these 
simple ideas. In particular, the absence of a MIT in 2D 
has been disputed by Finkelstein Q, who has empha- 
sized that corrections from interacting particles become 
important as T = is approached. These and more re- 
cent H considerations of interacting systems, have sug- 
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FIG. 2. Illustration that the effects of growth of 
L v oc 1/T P can be cut-off by a finite sample length L. 

gested that (3 can change sign for d — 2, i.e. that there 
can exist a 2D MIT. 

In fact, a wide variety of two-dimensional systems 
exhibit transitions like the MIT. Kosterlitz-Thouless- 
Berenshinskii transitions Q in 2D superfluids is an old 
example. There are more recent examples of proper 
conductor-insulator transitions in two-dimensional elec- 
trical transport problems. Granular films of supercon- 
ducting materials exhibit a transition from insulator to 
superconductor as the film thickness increases J?]]. The 
transition is driven by the superconducting fluctuations 
in competition with the Coulomb (charging) interactions 
between the grains of material [|| . In oxide superconduc- 
tors, the 2D sheets of oxygen atoms lead to a highly direc- 
tional order parameter and, depending on doping level, 
exhibit a superconductor to insulator transition some- 
what like that in granular films. (^] Between plateaux in 
the quantized Hall effect, which demands very high mag- 
netic field B such that the Landau energy hu> c ^> Ep, 
there is another 2D MIT manifest as a reversal of the 
temperature coefficient of the resistivity p. All 
of these examples appear in rather exotic circumstances 
(Landau quantization or Cooper pairing), but there is a 
manifestation of a 2D MIT in garden-variety conductors 
and that transition is the subject of the present work. 

In generic silicon metal-oxide-semiconductor field- 
effect transistors (MOSFETs), which comprise strictly 
2D electrons (or holes) between reservoirs of 3D carri- 
ers, Jl| there is a MIT H for MOSFETs with suffi- 
ciently low disorder and high carrier mobility [ p"5| . In the 
MOSFETs, the transition occurs among normal (non- 
superconducting) carriers at zero magnetic field. The 
transition occurs at relatively low electron densities n s 



of the order of 1 x 10 15 /m 2 , where the Coulomb energy 
(U ~ ^/n^) is larger than Ep . 

Clear signatures of the 2D MIT have been observed 
at B = in generic 2DES as a reversal of the sign of 
the temperature coefficient of conductance as the carrier 
density crosses through some critical value n c : dG(T) / dT 
changes from positive to negative. The same signature is 
found in a variety of materials including electrons in Si- 
MOSFETs, |L3],[l5) holes in GaAs/AlGaAs heterostruc- 
tures Jl7],[l8|], holes in SiGe/Si heterostructures ]lS| ], and 
electrons in Si/SiGe @. ' 

For all versions of the 2D MIT, a scaling of the resistiv- 
ity data against a control parameter can be accomplished 
in a form that directly follows from standard scaling con- 
siderations. In earlier work, this scaling form has been 
employed to explain the superconductor-insulator tran- 
sition in granular metal films Written in the form 
appropriate to describe the MOSFET experiment (where 
the control parameter is n s ), the scaling equation is ^ 

p{T,n s ) = f{\5 n \/T 1 ln = P{T/T Q ) 1 (2) 

where 8 n = (n s — n c ) /n c is a distance from the transition 
and z and v are exponents associated with the scaling of 
a microscopic correlation length £, which measures the 
regions over which the quantum interactions are effective 
- i. e. it is related to the phase coherence length L v . In 
the vicinity of the MIT, the correlation length is assumed 
to scale to infinity as £ = l/|c) rl | I/ . 

Note that standard scaling arguments are typically for- 
mulated in terms of the scale dependence of the conduc- 
tance at T = 0. It is clear that the T = limit of 
the theoretical models is a poor approximation to any 
experiment, since it presumes that the coherence length 
for the carrier excitations is infinite. As first emphasized 
by Thouless M), this assumption is false at any finite 
temperature; more recent arguments [B2| question its va- 
lidity even at T = 0. The predicted (3 can be related 
to more useful observables, by noting that the effective 
sample size for quantum interactions and interferences is 
not the patterned sample size but the length-scale L v 
on which carriers retain phase coherence or phase mem- 
ory. For most systems, this characteristic thermal length 
scales as a power of the temperature, viz. L v cx 1/T 1 / 2 , 
where z is the so-called dynamical exponent (z = 1/p 
in "weak localization" parlance). In this way, the scal- 
ing function in the length "domain" can be converted 
to T = -d{\ng)/d{\nT) in the T "domain" g. We 
emphasize that if the sample size L < L v , then temper- 
ature dependence may be "short-circuited" by the finite 
size effects as illustrated in figure H. 

The use of the quantum phase transition |2j| models 
to explain the 2D MIT carries a certain baggage with it. 
Such models are generally accepted for a superfluid (bo- 
son) ground state and transition only at T = to an in- 
sulating phase. For the quantized Hall effect or the gran- 
ular superconductor, there is an obvious choice of ground 



2 



state in the dissipationless (p — 0) behavior characteris- 
tic of these two physical systems. It begs the question, 
however, of the ground state in the high-mobility MOS- 
FETs: What is going on at T = ? The appearance 
of the 2D MIT in experiments has led to something of 
an avalanche of theoretical ideas, with "explanations" 
of the metallic behavior ranging from valley crossings, 
p4j triplet superconductivity |25| anyon superconduc- 
tivity |2(J , . . . There also have been suggestions that non- 
interacting models with spin-orbit scattering can explain 



the data |27|. These efforts nonwithstanding, the nature 
of the 2D metallic ground state and the 2D MIT remain 
some of the most challenging open problems of contem- 
porary condensed matter science. 



II. EXPERIMENT 

We have studied the conductance G of two different 
sets of high-mobility n-channel Si-MOSFETs. A set of 
large area (~ 1mm 2 ) FETs were fabricated on the (100) 
surface of silicon wafers doped with N a w 8.3 x 10 20 m -3 
acceptors. Corbino channels of length L = 0.4 mm and 
width w — 8 mm were formed with a poly-silicon gate 
above a 44 nm oxide layer. The measured residual ox- 
ide charge for these devices was w 3 x 10 14 /m 2 . These 
samples have been studied only at relatively high tem- 
peratures 1.2 < T < 4.2 K. Another set of samples made 
in a different process run on very similar starting mate- 
rial had a residual oxide charge of < 10 14 /m 2 . This run 
comprised samples of various lengths 1 /im< L < 256 /im 
and widths 11 fim< w < 500 /im. For the shorter sam- 
ples w > L, so that even though L is small enough 
to compete with important microscopic process length 
scales in the material there is some hope of inferring 
averaged properties as a result of having many "meso- 
scopic" elements in parallel. [H This later batch of 
samples have been studied at much lower temperatures 
0.01 < T < 4.2 K. All sets of FETs had rather high 
peak mobilities ~ 1 m 2 /V-sec, or they could achieve this 
regime with substrate bias. |ll| 

All measurements of G were conducted in an electri- 
cally shielded enclosure with standard lockin techniques. 
A source-drain voltage V s d was applied and the resulting 
current I s d was recorded as a function of the temper- 
ature T, the magnetic field B, or the gate voltage V g , 
which controls the carrier density n s = C ox (V g — Vth)- 
Vth is the threshold for populating the channel, which can 
be inferred from higher temperature transconductance or 
from Shubnikov-de Haas measurements at low tempera- 
tures, and C ox is the specific capacity of the gate oxide. 
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FIG. 3. (a) a(T,n s ) of a large MOSFET for T = 1.32 
(o), 1.84 (+), 2.30 (v), 2.70 (□), 3.00 (A) and 3.50 K (O). 
(b) The same data re-scaled according to Eq. ^ to form 
a single function with two branches (bottom branch for the 
insulating phase and top branch for the metallic behavior) (c) 
The data re-scaled again in accord with Eq. Q so that the 
symmetry for insulating and metallic conductance is evident 
up to T = To, whch is marked by the dashed line. The dotted 
line marks the conductivity a c = 0.8e 2 /h at the MIT. (These 
data were recorded at finite magnetic field B = 0.5 T, but 
this has no bearing on the shapes of any of the curves for this 
sample.) 
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III. RESULTS FOR B — 

Our experiments demonstrate first that the disorder 
in the devices must be low relative to other important 
energy scales. For a moderately disordered sample, we 
find that the slope of G(T) is consistently positive: all 
values of n s scale as insulators and the conductivity van- 
ishes as T — > 0. For a substantial back-gate bias {e.g. 
V su b — — 9 V for the sample shown in Fig. ||), how- 
ever, the mobility increases enough to permit a MIT at 
a value of n s — n c ~ 1.65 x 10 15 /m 2 . The MIT appears 
as a change of sign of the slope of cr(T) = (L/w)G(T) 
at n c (see Fig. 0(a)). We infer that the peak mobility 
has increased by about 20%, but this clearly is sufficient 
to change the behavior of the 2DES dramatically, even 
though the value of a at n s — n c has not changed much. 

More importantly, all of the conductance curves 
G(V g ,T) = G(n s ,T) can be scaled (see figure |(b)) 
against 5 n — (n s — n c )/n c to form a single two-branched 
function [Fig. ||(b)] as expected from the scaling argu- 
ments mentioned above. We take this as a signature of a 
quantum phase transition in the 2DES. In figure |3](c) we 
plot 



a(S n , T) = a c exp (^A5 n /T 



l/z 



(3) 



for all data from figure §(b), and find agreement with 
the theoretical prediction |5j for the temperature depen- 
dence in the quantum critical region T > To, where the 
crossover temperature To oc |5 n | ziy is shown by the dashed 
line in Fig. ||(c). 

Very similar scaling of a has been obtained in other 
2DES samples in a variety of physical circumstances. 
lljU^HI! A related scaling of a(E,T,n s ), where E 
is the electric field applied between the source and drain, 
has been observed in some 2DES. (li) This collection of 
experimental results is overwhelming evidence that the 
MIT is a quantum phase transition, and exhibits simi- 
larities with transitions in granular superconducting thin 
films, oxide superconductors, and plateaux transitions in 
the quantized Hall effect. 

In shorter MOSFETs we find a similar scaling of a for 
moderate temperatures 1 < T < 4.2 K. Our results are 
illustrated in figure |] which contains measurements from 
a L = 1.25 fim, w = 11.5 MOSFET. The reversal of 
sign of the slope dG/dT occurs at essentially the same 
value as for the larger MOSFETs, and the conductivity 
<7 C at the MIT is ~ e 2 /2h in agreement with many other 
experiments on MOSFETs. A similar value of a c ob- 
tains in completely different physical circumstances such 
as the granular films or the quantized Hall systems, but 
there the amplitude of a is governed by physics (charging 
energy and phase fluctuations or the dominance of a par- 
ticular momentum mode) that are not directly germane 
to the MIT. 
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FIG. 4. Scaling of a(n s ,T) for a 1.25 /im long 11.5 /im 
wide MOSFET for T = 1.97 (A), 2.48 (o), 2.68 (□), 3.06 
(x), 3.53 (v) and 3.93 K (o, +). The dashed line marks 
a c = 0A5e 2 /h. 
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FIG. 5. Scaled a(n s ,T) for T = 4.2,1.2 and 0.4 K; 
Vsub = 10 V; n c « 1.8 x 10 15 /m 2 . The dashed line marks 
o c = 80e 2 /h, and the dotted lines guide the eye. 



Yet a third set of samples with electrons confined at the 
interface between Si and SiGe has exhibited very similar 
20 1 This result (see figure ||) is especially in- 



response. 

triguing, because the mobility is very high compared with 
the MOSFETs experiments |^9| - in the same range as for 
electrons in GaAs heterostructures, which have not ex- 
hibited the MIT to date. Again the slope dG/dT changes 
sign at a concentration n c = 1.8xl0 15 /m 2 , but owing to 
the higher mobility, the value of a c ~ 80e 2 /h is very 
large compared to the MOSFET experiments. This re- 
sult proves the existence of the MIT in yet another ma- 
terials system. More interestingly, it lays to rest once 
and forever the suggestion (quoted widely and believed 
even more widely) that the value a c is a "universal" num- 
ber. Indeed, subsequent studies of the MIT in a 2D hole 
systems in SiGe/Si [[t9| and GaAs/AlGaAs heterostruc- 
tures Jl8| have found a c ~ 2e 2 /h, almost a factor of seven 
larger than in some Si MOSFETs. 

A much more plausible expectation is that r s = U / Ep 
will be a universal number. This, in fact, is not born 
out across all samples: r s ~ 12 — 19 at the transi- 
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tion for electrons in Si MOSFETs and Si/SiGe quantum 
wells, r s ~ 13 for holes in SiGe/Si heterostructures and 
r s ~ 11 — 23 for holes in GaAs/AlAs heterostructures. 
These numbers, however, are much closer to each other 
(within a factor of two) than the values of ct c , which are 
distributed across more than two orders of magnitude. 

In general, one expects the scaling exponents, such as 
z and v, to reflect only the symmetry of the problem in 
question, and thus to be universal numbers. In the ex- 
periments, however, a large range of exponents are found. 
In the Si-MOSFETs at B = 0, zv = 1.6 ± 0.2 has ap- 
peared in two different sets of samples. |13p5f| The same 



exponent has been reported for holes in SiGe/Si quan- 
tum wells. [H For holes in GaAs/AlAs heterostructures 
flisf or electrons in Si/SiGe heterostructures |2(J, how- 
ever, zv > 4. It is worth noting that these two experi- 
ments were on samples with very low disorder (hp I 3> 1) 
by comparison with the MOSFET and p-SiGe samples 
where the relative effect of disorder is much stronger 
(k F l ~ 1). In the very short MOSFETs |oH|, zv is 
even larger - zv = 16 ±4 for L w 1/im. Furthermore, the 
range of "scalable" temperature dependence is smaller 
than in the larger samples. In fact for \8 n \ > 0.15 on the 
insulating side, the lowest temperature curve (A) is al- 
ready failing to scale with the higher temperature data. 
In the sample of figure [|, the data obeyed Eq. (||) well 
for 1.2 < T < 3.6 K (T = 1.2 K was the lowest tem- 
perature studied in that experiment), and for other large 
MOSFETs [jl3| the range of temperature extended down 
to 0.05 K. We attribute this saturation of the T depen- 
dence to the cut-off of L v and hence £ by the sample 
length L. 



IV. MAGNETOCONDUCTANCE 

An applied magnetic field B alters the behavior of the 
carriers in the 2DES in different ways depending on the 
angle between the field and the plane of the 2DES. The 
perpendicular component of the magnetic field tends to 
bend the path of the carrier (the Lorentz force tends to 
make the carriers form circular orbits), as well as induc- 
ing a splitting of the energies of different spin orienta- 
tions and consequently to polarize the spins of the carri- 
ers out of the plane of the 2DES. A parallel component 
of the magnetic field probably has little effect on the or- 
bital motion of the carriers, but it still splits spin ener- 
gies and it polarizes the spins into the plane. Studies of 
G(B) have been employed fruitfully for decades as probes 
of the transport mechanisms of disordered systems. |3S| 
There are detailed calculations of the functional form of 
G(B) for various conditions among non-interacting car- 
riers p^ |. The contributions from Zeeman interactions 
(triplet terms in Hartree approximation) have been cal- 
culated and for \igB < ksT (pg is the carrier magnetic 
moment and ks is Boltzmann's constant) lead to a neg- 
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FIG. 6. Magnetoconductance for a large MOSFET with 
Vsub = -9 V and n c = 1.65 x 10 15 /m 2 . The different curves 
arise from the values of n s written on the plot (units of 
10 15 /rn 2 ). 



ative parabolic magnetoconductance in a perpendicular 
magnetic field. In contrast, quantum interference of car- 
riers (orbital effects) can lead to positive slopes of G(B) 
near B — 0. 

The MIT in MOSFETs is quenched by a large paral- 
lel magnetic field (B > 1 T) |3^,Q, which in turn im- 
plies that polarizing the electron spins in the plane of the 
carrier motion has a dramatic effect on the correlations 
of the carriers. The conductance of nominally metallic 
2DES (i.e. n s > n c ) decreases by orders of magnitude 
and the slope of G(T) reverses from negative to positive, 
indicative of insulating samples. The application of B 
perpendicular to the plane of the 2DES also quenches the 
metallic phase |3^,|3(|. Moreover, careful measurements 
of magnetoconductance in the quantum critical region in 
the presence of a perpendicular B provide quite a bit of 
insight into the relative importance of spin interactions 
and orbital motion at the MIT. (IB) 

Figure contains representative plots of magnetocon- 
ductance (MC) Aa/a(B = 0), where Act = a(B) - 
a[B = 0). Clearly, there is a positive magnetoconduc- 
tance and a negative magnetoconductance contribution 
to each curve, and the positive contribution is more ac- 
centuated near n s = n c = 1-65 x 10 15 /to 2 . As we will 
see below, this is because the negative contribution has 
dropped to a minimum value at this point. 

Each curve in figure o can be written 



Aa/a(0) = pf(B) - qB 2 



(4) 



where p, q > 0. The assumption of parabolic form for 
the negative term is justified because (1) it fits the data 
and (2) it is the form expected for the electron-electron 
interaction contribution for B < ksT/fig ~ 0.8 T [ ^7| . 
We have decomposed the data following formula (^)and 
find that the positive component pf(B) of the magneto- 
conductance does not depend on n s within the scatter of 
data. These contributions are plotted in figure ^(a) and 
fall atop one another for B < 0.8 T. The coefficient p is 



5 



0.12 - 




0.0 0.2 0.4 0.6 0.8 1.0 



B (T) 



0.10 



I 

(b) 




B = 1 T / 

/ 

/ 

/ 

/ 

/ 

/ — 

. / 
• • 


- • 




-O 

I 


0~ 


— --o-°-— °- 
o 



1.35 1.45 1.55 1.65 1.75 1.85 
n s ( 10 15 m " 2 ) 



FIG. 7. (a) The positive contribution pf{B) to 

Aa/a(B = 0) for all 1.40 < n s < 1.80 x 10 15 /m 2 . (b) 
Amplitudes p (o) and q (•) of Aa/a(B = 0) as a function n s 
for B = 1 T. The lines guide the eye. 

therefore a constant. The coefficient q, i.e. the electron- 
electron interaction contribution to MC, exhibits a min- 
imum near n s — n c as is obvious in figure @(b) . We have 
found exactly the same behavior of MC in short MOS- 
FETs |30| , |38|] at temperatures down to 0.04 K, in spite of 
the presence of sizable conductance fluctuations. 

Both positive and negative MC have been observed in 
other MOSFETs as well @. Since orbital effects (posi- 
tive magnetoconductance) are absent in an applied par- 
allel field, our results show clearly that MC of a 2DES 
does depend on the angle between the magnetic field and 
the plane of a 2DES, at least in the quantum critical re- 
gion. Away from the critical region, the orbital effects 
at low fields become small compared to spin effects [see 
Fig. 0(b)], and only negative MC is observed. Therefore, 
a recent claim |36[ ] about the absence of any angular de- 
pendence of MC is valid only far from the critical regime, 
in the metallic phase where the experiment p6[ was car- 
ried out. 

V. SCALING OF THE CONDUCTANCE 

We can extract (3t oc (3 from our data by using the 
assumption of power-law dependence of the correlation 
scale (the phase coherence length) on T. |?2|,|| One 
such plot for four of our samples is given in figure ||. 
Anticipating the dependence of L v suggested in the in- 
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FIG. 8. Experimental values of f3r scaled by l/zu for four 
samples: two large area (O an d □) MOSFETs and two short 
(V and O) samples (for high enough T). 

traduction, we have scaled the bare (3t multiplying it by 
zv. The two large MOSFETs have zv = 1.6, and the 
two short MOSFETs yielded zv = 16. In spite of the 
disparity, all four curves produce the same scaled scal- 
ing function indicating that the underlying (3 is the same 
for all samples. A noteworthy feature is that (3 is essen- 
tially linear as it crosses through j3 — as expected from 
recent scaling arguments jH). For non- interacting parti- 
cles, (3 < as illustrated by the dotted lines in figure 
|l|. Our experiment (which corroborates other measure- 
ments |^]) proves that (3 extends into the metallic region 
and thus provides clear evidence as to the inadequacy of 
noninteracting models of charge carriers. 

An even more striking feature is evident from the 
present data: [3 decreases towards zero at large <?, as 
proposed in Ref. M. Although this is by no means a 
necessary condition for such an exotic metallic state, it 
seems as if the system tries to restore Ohm's law in the 
large conductance (weak disorder) limit. 

VI. SUMMARY 

Our experiments corroborate other experiments on 
2DES and prove the existence of a metal-insulator transi- 
tion in two-dimensions for sufficiently low disorder. The 
conductances can be reduced to a single (two-branched) 
scaling function that is consistent with recent proposals 
for Coulomb driven quantum phase transitions. Mea- 
surements for different sample lengths and different lev- 
els of disorder demonstrate certain universal features of 
the transition. The transition occurs near a certain value 
of n c for a given density of states (carrier effective mass). 
The conductance at the transition is not a universal con- 
stant. Magnetoconductance experiments suggest that 
spin polarization and Coulomb interactions of the car- 
riers have a dramatic effect on the transition. The prod- 
uct of conductance scaling function (3t and the correla- 
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tion scale exponents zv shows that the underlying scaling 
function f3 is the same for all MOSFETs. The inferred j3 
violates the predictions from non-interacting models and 
supports recent predictions based on interacting carriers. 
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